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Abstract 

Periodicity and relaxation are investigated for the trajectories of the states in 
one-dimensional finite cellular automata with rule-90 and 150. The time evo- 
lutions are described with matrices. Eigenvalue analysis is applied to clarify 
the maximum value of period and relaxation. 
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I. INTRODUCTION 



Cellular Automata are one of the simplest mathematical models for nonlinear dynamics 
to produce complicated patterns of behaviour, which had been originally introduced by 
von Neumann [jT|. Wolfram had reintroduced cellular automata as a model to investigate 
complexity and randomness |0. He investigated many fundamental features of them 
Since then many authors have made efforts to clarify the properties of cellular automata 
and applied to natural systems p. 

One-dimensional cellular automata are described by the discrete time evolution of site 

af. 

ai{t + 1) = F[ai-r{t),ai-r+i{t), ai{t), ai+r(t)], (1.1) 

where takes k discrete values over Z^. The simplest model, elementary cellular automaton, 
consists of sites with two internal states over Z2 interacting with the nearest neighbour sites 
(r = 1). Wolfram introduced a naming scheme for these models and classified the behaviour 
of cellular automata into four classes [1,0. 

Most authors have worked on cellular automata within the scope of the infinite number 
of sites. A few works have concerned periodic boundary condition (cylindrical automata) 
P|7HT0|. In our previous paper |Tl]] (referred as paper I) we had investigated the periodic 



orbits of finite rule-90 cellular automata with Dirichlet boundary condition. We analyzed 
the eigenvalue equations of the transfer matrices which describe the time evolution of the 
system. In the present paper the method is extended to the rule- 150 case. Some results 
obtained in paper I will be cited again for completeness. 



II. MODELS AND NUMERICAL RESULTS 

We review the numerical results obtained in paper I on so-called rule-90 cellular automata 
following Wolfram's naming scheme. The time evolution of the i-the site ai{t) G {0, 1} 
(i = 1 ~ A^) is described as a sum modulo 2 of the nearest neighbour sites: 
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ai{t + 1) = ai_i(t) + a,+i(t) mod 2. (2.1) 

We use the Dirichlet boundary conditions, = a^+i = 0. In the Wolfram's classification 
the rule-90 model belongs to the third class which shows the chaotic behaviour. The time 
evolution is also expressed by matrix: 

A{t + l) = UA{t), (2.2) 

where A{t) = *(ai(t), a2(t), . . . , aAr(t)) describes the state at t and the transfer matrix U is 
given by 

rl j=^±l, 

U., = (2.3) 
otherwise. 

In paper 1 we had found the periodic structures of the transfer matrices U numerically as 
shown in Table |I[ They are summarized as follows: 

U^^=I, (iV is even), (2.4a) 
IfU^+TTM = jj^N^ (jv is odd), (2.4b) 
= 0, (iV = 2" - 1). (2.4c) 

For cellular automata with even number of sites (Eq. ( |2.4a] )), every state except the null 
one (all sites are zero) is on the orbits with period not exceeding Hat. The concrete periods 
depend on the initial states. The period Hat corresponds to the least common multiple of 
them, and we call it the maximum period. The states on the orbits with the maximum 
period have the lowest symmetry. The states with some symmetries are on the orbits with 
shorter periods than n^r. For cellular automata with odd number of sites except = 2" — 1 
(Eq. (|2.4b|) ), some states belong to the orbits with period Hn or less and the others except the 
null state are drawn to the periodic orbits after some time steps not exceeding the maximum 
relaxation ttat. It is very interesting that in the N = 2^ — 1 {n E Z) case every state is drawn 
into the null state after at most steps (Eq. ( p.4c| )). In this case the configuration space 



has only one basin with the null state at the center. Schematic features of trajectories are 
shown in Fig |l[ 



Next we investigate the rule- 150 cellular automata, which also belongs to the third class 
in the Wolfram's classification. The time evolution is described as a sum modulo 2 of the 
site itself and the nearest neighbors: 

ai(t + 1) = aj_i(t) + ai(t) + aj+i(t) mod 2. (2.5) 

The transfer matrix is given by 

{1 j = i±loTj = i, 
(2.6) 
otherwise. 

The periodicity of the transfer matrices are also found numerically as shown in Table |. There 
are no cases drawn into the null state as described by Eq. ( p.4c| ) . For cellular automata with 
N = 3n + 2 there appear the periodic orbits of period not exceeding IItv with the relaxation 
path whose maximum length is 71^. For the others 7^ 3n + 2 all states are on the periodic 
orbits except the null state. 

The maximum periods IIjv of the rule-90 and rule- 150 cellular automata are found to 
coincide each other in many cases of the number of sites (See Table |, Figs. ^ and ^. 
Grassberger had reported the similarity between the rule-90 and rule-150 cases in behaviour 
T^. Especially the steepest peaks N = 6,10,18,22,28 coincide between the rule-90 and 



rule-150 cases, whose periods are expressed by Hat = 2^/^"*"^ — 2. The maximum periods 
grow exponentially for the average of all number of sites except = 2" — 1 cases, in which 
all state are drawn into the null states within steps for rule-90 and n^r = + 1 for 
rule-150. 



III. EIGENVALUE ANALYSIS 

The trajectories of rule-90 cylindrical automata had been investigated by Martin, 
Odlyzko and Wolfram They analyzed characteristic polynomials which describe the 
states of cellular automata. The method seems to work well in the periodic boundary cases. 
To investigate the periodic structures analytically in the Dirichlet boundary cases, we had 
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introduced the simpler method to analyze the eigenvalue polynomials in paper I. The eigen- 
value equation of the transfer matrix 

UA = -XA, (3.1) 

reads the secular equation: 

D^{X) = \U + XI\ = 0. (3.2) 

All calculation are carried out over Z2 since a site takes binary values. 

The recursion relation of the eigenvalue polynomials for rule-90 cellular automata 

D^(A) = XD^'-^X) - D'^-^X), (3.3) 

gives the explicit form of D^{X). The eigenvalue polynomial D^{X) is over Galois field of 
order 2, GF{2), namely all coefficients of A* are over Z2: 

J max 

D^(A) = £ (Cf mod 2) A^-2^ (3.4) 

j=0 

where jmax = [N/ 2J , [J is a Gaussian symbol and 

Cf^{-iy(^^-^y (3.5) 

Note that the definition of is slightly different from that in paper I. 

The eigenvalue polynomials are given by the replacement A — > A + 1 in Eqs. (|3.3| ) and 
I4D for the rule-150 cellular automata: 

D^(A) = (A + l)D^-i(A) - D^'-^X), (3.6) 



Z}^(A) = J2 (Cf mod 2) (A + lf-^\ 

N l{N-k)/2\ 

= E E (ciJmod2)A\ (3.7) 

(3.8) 



Some examples of D^{X) for the rule-90 and rule-150 cellular automata are shown in Tables H 
and g. 

These eigenvalue equations enable us to find the maximum periods and maximum length 
of relaxation path to the periodic orbits. First we study nilpotent cases, 
paper I we found that those happen for the rule-90 with = 2" — 1 sites and P{N) = N. 
All states are drawn into the null state within steps or less. The rule-150 model does not 
show the similar behaviour. 

Next we discuss the cases that there are constant terms in polynomials. Those happen 
for the rule-90 with even number of sites since = (— l)^/^"^ f^j, even A^. For the rule- 

150 case the polynomials has a constant term for A^ 7^ Sri + 2 by Eq. (|3.6|) . The eigenvalue 
polynomials reduce to the simple form as A^*-^-* + 1 = by multiplying some power of A 
and repeatedly substituting the eigenvalue equation to itself ||T3| , p!^ . The minimum value 
of the power P{N) corresponds to the maximum period Un. If the eigenvalue polynomial 
is factorized, we are able to get shorter periods depending on the initial state from those 
factors by the procedure mentioned above. For instance, we consider A^ = 4 rule-150 cellular 
automata (Fig. ^). The eigenvalue equation A^ + A^ + 1 = reduces to A^ + 1 = and the 
maximum period is 114 = 6. The eigenvalue equation is also factorized to (A^ + A + 1)^ = 0. 
We find another solution A^ = 1 from A^ + A + 1 = 0. Therefore there are period 6 and 
period 3 orbits for A^ = 4 rule-150 cellular automata. 

Explicit expression of the maximum period II^v is simply obtained for A^ = 2" — 1 rule-150 
cellular automata. In this case the identity (^~"') mod 2 = holds for j 7^ 0. The non-zero 
elements of C'JJ. are Cq^ = (^^^ . By the identity mod 2 = 1 for all non-negative 

integers k, the eigenvalue polynomial reduces to D^(\) = Y.k=o^''- Following the above 
procedure, we find the maximum period IItv = A^ + 1- 

Finally we study the cases that the lowest powers of polynomials are greater than 0, 
namely they have the forms as D^{\) = A*''-^-* ■ /(A), where /(A) is a polynomial with a 
constant term. The number of null solutions, p{N), corresponds to the maximum length of 
relaxation path ttat. Applying the above procedure to /(A), we obtain the maximum period 
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and shorter ones. For example, we show the case of = 5 rule-90 cellular automata (Fig. 
The eigenvalue polynomials is D^{X) = A ■ (A^ + 1). The maximum length of relaxation 
path is TTs = 1. From A^ + 1 = 0, we find the maximum period lis = 4. We also find shorter 
ones 2 and 1 by factorizing A'' + 1 = (A^ + 1)^ = (A + 1)^. 

A simple recursion relation of the maximum period and maximum length of relaxation 
are found for rule-90 cellular automata with odd number of sites. The new recursion relation 

D^(A) = X^D^'-^X) - D''-\X), (3.9) 

holds by virtue of algebra in GF{2). This gives a simple form D^"+^(A) = Ai5"(A^). The 
maximum period and the maximum length of relaxation path of = 2n + 1 rule-90 cellular 
automata are described by those of A^ = n sites as 

n2„+i = 2n„ + i, (3.10) 

TTan+i = 27r„ + 1. (3.11) 



IV. CONCLUDING REMARKS 

We investigate the periodic structure of one- dimensional rule-90 and rule-150 cellular 
automata with Dirichlet boundary condition. We find three types of behaviour. The first is 
periodic one which appears in cellular automata with even number of sites for rule-90 and 
A^ 7^ 3n + 2 for rule-150. The second appears in the cases of odd number of sites for rule-90 
and A^ = 3n + 1 for rule-150. There are some periodic orbits and irreversible relaxation 
paths to them. The peculiar behaviour happens to the case A^ = 2" — 1 for rule-90. All 
states are drawn into the null state within A^ steps. 

The eigenvalue equations of the transfer matrices are analyzed. The maximum period is 
obtained by finding the minimum power to satisfy A^^^-* -|- 1 = 0. Shorter periods are given 
by factorizing the eigenvalue polynomials. The number of null solutions of the polynomials 
gives the maximum length of relaxation path. For some special cases, we find the explicit 
forms of the maximum period. Distribution of periods and relaxation is still not clear. 
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Roots of the eigenvalue equations, in general, are found over the Galois extension 
of finite field p^JT^ . For all positive integer N there are primitive polynomials V{x) 
of degree N over GF{2). One of the roots, a, of V{x) generates the Galois exten- 
sion GF{2^), whose elements are {0, 1, a, a^, . . . , and a"^^^^ = 1. Other roots of 
V{x) are called conjugate of a and generate the isomorphic Galois extensions. A gen- 
eral polynomial of degree has roots over GF{2^). A eigenvalue polynomial has A^ 
roots {ttj} G GF{2^) {i = 1,2, . . . , N). For instance A^ = 4 rule-150 case, the equation 
+ + 1 = has roots in {0, 1, a, o? , . . . , a^^} (a^^ = 1), where a is one of roots of forth 
primitive polynomial A^ + A + 1 = 0. Explicitly the roots are oc' and o^^ with multiplicity 2 
for each. The minimum common multiplier of the roots, which satisfies = 1, is a^^ . This 
seems to suggest the maximum period 6, namely (a^)^ = [ot^^Y' — '^'^ — 1- This procedure, 
however, does not work for other cases. For example, the eigenvalue polynomial of A^ = 3 
rule-150 case is A^ + A^ + A + 1. It can be f actor ized to (A + 1)^ and the root is A = 1 with 
multiplicity 3. The periods, however, are 4, 2 and 1. More number theoretical studies are 
expected. 

The elementary cellular automata are also a subject to build a built-in self-test of VLSI 
T^,|TB|. Usually the shift-registers which generate pseudo-random sequences of length 2^ — 1 



with A^ registers, are used for build-in self-tests. Our results show the elementary cellular 
automata can also generate exponentially long but not maximum sequences. By fine-tuned 
mixture of rule-90 and rule-150 cases, hybrid cellular automata, it is shown to be able to 
produce a maximum length pseudo-random sequences. 
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FIGURES 

FIG. 1. Schematic features of the trajectories of cellular automata: (a) simple periodic orbit, 
(b) periodic orbit with relaxation, (c) limit point. 

FIG. 2. Distribution of the maximum periods of the rule-90 cellular automata. The bold 
line is the average of the periods except N = 2^ — 1 cases. The broken one denotes the curve 
Un = 2^/2+1 - 2, which fits peaks of iV = 6, 10, 18, 22, 28. 

FIG. 3. Distribution of the maximum periods of the rule-150 cellular automata. The bold line 
is the average of the periods of all number of sites. The broken one denotes the same as in Fig. ^. 

FIG. 4. The orbits of = 4 rule-150 cellular automaton. All states are classified into 3 orbits 
excepts the null states. Two of them are period 6 and the other period 3. The null state is an 
isolated fixed point. The periods is same as those for rule-90, though the detail behaviors are not 
(See Fig. 1 in paper I). 

FIG. 5. The orbits of = 5 rule-90 cellular automaton. This figure is the same as Fig. 2 in 
Paper I. 
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TABLES 



TABLE I. Periodicity of the transfer matrix U for rule-90 and rule-150 cellular automata. 

rule-90 rule-150 

f/3 = ^^ = 1 

U^ = I U^ = I 

= u = u^ 

U^ = 7 U^ = I 

= / U^^ = C/2 

= U U^^ = 1 

i762 = 1 U^'^ = I 

C/ii = C/3 = C/8 

f/126 = / U^^ = I 

= U U^^ = I 

C/30 = / U^^ = C/2 

= 7 J7^6 = I 

i730 = 7 = 7 

C/29 = [/ [732 = C/4 

[yl022 _ J jyl022 _ j 

U^^ = = I 

= / [/128 ^ f/2 

^4094 _ J ^4094 _ j 

C/23 = C/7 ^724 ^ [/16 

fj2046 _ J jy2046 _ j 

^1022 _ J ^1024 _ ^2 
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27 U^^ = C/3 [7^6 ^ J 

28 [/32766 _ J ^32766 _ j 

29 C/^^ = [/ [/^^ = C/^ 

30 ^ J ?762 ^ / 

31 = U^'^ = 1 

32 = 7 = C/2 
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TABLE II. Eigenvalue polynomials for rule-90 cellular automata. 



A3 

+ A2 + 1 
A^ + A 
A^ + A^ + 1 
A^ 

A^ + A^ + A^ + 1 

A^ + A^ + A 

A^o + A8 + A^ + A2 + 1 

All + 

Ai2 + A^o + a8 + A2 + 1 

Ai3 + A9 + A 

Ai4 + + a8 + 1 

Ai5 

Ai6 + Ai4 + Ai2 + a8 + 1 
A^^ + Ai3 + a9 + A 

Ai8 + Ai6 + Ai2 + A^o + AS + A2 + 1 

Al9 + All + 

A^o + Ai8 + Ai6 + Alo + A8 + A^ + a2 + 1 
A21 + Ai7 + A9 + A^ + A 

^22 + ;^20 + Al6 + A8 + A6 + A^ + 1 

A23 + A^ 

A24 + A22 + A20 + Ai6 + A6 + A^ + 1 

A25 + A^i + Ai^ + A5 + A 

A26 + A24 + A20 + Ai8 + Ai6 + A^ + A2 + 1 

A27 + Ai9 + A3 
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28 + + + + Ai6 + a2 + 1 

29 A29 + A25 + Ai7 + A 

30 A^o + A28 + A24 + A^e + 1 

31 A^i 

32 A32 + A30 + A28 + A24 + Ai6 + 1 
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TABLE III. Eigenvalue polynomials for rule-150 cellular automata. 



+ A2 + A + 1 
A^ + A2 + 1 
A5 + A4 
A^ + A^ + 1 

A^ + A^ + A^ + A^ + A^ + A2 + A + 1 

A8 + A6 + A^ 

A^ + A^ + A^ + A^ + A + 1 

A^o + A^ + 1 

All ^ ;^10 + ^9 + A8 

Ai2 + A^o + A8 + A^ + 1 

Ai3 + + + A^ + A + 1 

Ai^ + Aio + A8 + A6 + A2 

Al5 + Al4 + Al3 + Al2 + All + ;^10 + + + 

+A6 + A^ + A^ + A^ + A^ + A + 1 
Ai6 + Ai4 + Aio + A8 + A^ + A2 + 1 
Ai7 + Ai6 + Ai3 + Ai2 + A^ + A^ 
Ai8 + Ai2 + Aio + A8 + A^ + A2 + 1 

Al9 + Al8 + Al^ + Al6 + All + ;^10 + 

+A8 + A^ + A^ + A + 1 
A^o + Ai8 + Ai6 + Aio + A^ 
A^i + A20 + a9 + A8 + A + 1 
A22 + Ai8 + Ai6 + AS + 1 
A23 + A22 + A^i + A20 + Ai9 + Ai8 + Ai7 + Ai6 

A24 + A22 + Al8 + A8 + 1 

A25 + A24 + A21 + A20 + Ai7 + Ai6 + A^ + A8 + A + 1 
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26 + + + A^o + A^ 

27 A27 + A26 + A25 + A24 + All ^ ;^10 + 

+A^ + A3 + A^ + A + 1 

28 A28 + A26 + A24 + A20 + Ai6 + Ai2 + Aio + A^ 

+A4 + A2 + 1 

29 A29 + A28 + A^i + A20 + Ai7 + A^^ + A^^ + A^^ + A^ + A^ 

30 A^o + A26 + A24 + a22 + Ai8 + Ai4 + A^^ + A^ 

+A6 + A^ + 1 

31 A^l + A^O + A29 + A28 + A27 + A26 + A25 + A24 + A23 

+ A^i + A^o + Ai9 + Ai« + Ai^ + Ai6 + A^^ 

+A14 + Al3 + Al2 + All + ^10 + ^9 + ^8 + ^7 

+A6 + A^ + A^ + A^ + A2 + A + 1 

32 A32 + A^o + A26 + a24 + A22 + Ai8 + Ai4 + Aio 

+A8 + A6 + A2 
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Fig.2 S. Tadaki Orbits in one-dimensional finite linear cellular automata 
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Fig.3 S. Tadaki Orbits in one-dimensional finite linear cellular automata 




Fig.4 S. Tadaki Orbits in one-dimensional finite linear cellular automata 
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Fig.5 S. Tadaki Orbits in one-dimensional finite linear cellular automata 



